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Summary 

A review of turbine-blade coolant-side heat-transfer 
technology is needed as inaccuracies in many of the 
correlations or extrapolations of correlations can result in 
large errors in the prediction of turbine-blade metal tem- 
peratures and life. Designers have often used correlating 
heat-transfer and flow equations because of a general 
acceptance of its usage and do not have the time to 
inquire about the inaccuracies of the correlations, the 
conditions, and the range of applicability. 

This review indicated that even for ideal conditions the 
correlations for heat transfer and pressure drop were at 
best about ±10 percent accurate and some approached 
±35 percent error. Further research is needed to provide 
needed data and improved accuracy of predictions for 
local turbulent heat transfer and pressure drops for inlet 
regions of passages, noncircular passages with and with- 
out turbulators (ribs and pins), impingement cooling 
coupled with film injection of the coolant, and, particu- 
larly, for passages under conditions of rotation. 


Introduction 

Much information has been obtained on the heat- 
transfer coefficients and pressure losses of airflow 
through passages. More continues to be obtained to 
extend the range of conditions, determine the effect of 
variations of cooling methods and configurations, and 
improve the accuracy of the equations that correlate the 
data. Those working directly with cooling of turbines 
have organized and condensed much of the information 
and have selected correlations they consider most appro- 
priate to their usage and cooling methods. However, 
much of this information is proprietary. Those who do 
not have access to the proprietary data must rely on 
published literature. Designers generally use correlating 
equations because of the general acceptance of their 
usage, and often extrapolate these equations beyond their 
range of applicability. A review of coolant-side heat- 
transfer technology is needed as inaccuracies in some of 
the correlations used or extrapolations of these correla- 
tions can result in large errors in the prediction of 
turbine-blade metal temperatures and life. Even a ±10 
percent inaccuracy in the coolant-side heat-transfer co- 
efficient can mean an uncertainty of ± 15 K (27° F) in the 
surface metal temperature of an advanced turbine-blade 
design (ref. 1). This uncertainty in metal temperature can 
affect blade life by a factor of about 2. 

The purposes of this study were to (1) review some 
published literature on heat transfer and pressure losses 
of air flow through passages for several cooling methods 
generally applicable to gas turbine blades, (2) select the 


more useful correlating equations for coolant passages in 
nonrotating and rotating blades, (3) assess the status of 
turbine-blade internal air-cooling technology, and (4) 
describe areas where further research is needed. The 
cooling methods considered include convection cooling, 
impingement cooling at the leading edge and at the mid- 
chord, and convection cooling augmented by pins and 
internal ribs. Also in this study, emphasis was placed on 
turbulent flows because turbulent flow is generally 
needed to promote the high degree of blade cooling 
required in gas-turbine engines. Transpiration cooling 
and film cooling, although important to the cooling of 
advanced turbines, were not considered in this study. 
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flow area in pin array section, = v/L 
ratio of jet hole area to opposing impingement 
surface area (open area ratio) 
correlation function for impingement heat 
transfer on flat plate (eq. (37)) 
correlation constants (eqs. (34) to (36); also see 
table IV) 


dimensionless quantity, B = 


(V2 A* 0 C d L)/z 


correlation function for impingement heat 
transfer on flat plate (eq. (37)) 


width of an equivalent two-dimensional nozzle 
having equal flow area, = xe?2/4c; correlation 
constant (eqs. (34) and (35); also see table IV) 
correlation constant 


coefficient of discharge 

center-to-center distance between impingement 
holes; correlation constant (eqs. (34) and 
(35)); also see table IV) 
specific heat 
pipe inside diameter 

hydraulic diameter =4A/(P; for impingement 
cooling in leading edge cavity, 

= 2 bc(n - 1 )/[b + c(n - 1)] 
characteristic diameter as defined after eqs. (51) 
and (54) 

effective pipe diameter =4//7 t 
impingement hole diameter; pin fin diameter 
height of turbulators 
roughness Re ynol ds number, =eu’/ v 

= (e>AD)ReV/72 

friction factor 

mass flow rate per unit area, = W/A 
heat-transfer funct ion , 

= I//(2St)-l]/V7/2 + K e + 
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g =g Pr-0.57 

H enthalpy 

h heat -transfer coefficient 

ij correlation exponents (eq. (47)) 

K geometry factor for turbulent flow; total pres- 

sure loss coefficient 
k thermal conductivity 

L length of flow passage; length of pin array 

section or jet -impingement plate in streamwise 
direction 

l half arc length of leading-edge cavity from 

stagnation point; pin fin height 
M Mach number; initial crossflow to total jet flow 

ratio, =W C /Wj 

m correlation function or constant for eqs. (33), 

(37), and (48) 

Nu Nusselt number, = hD/k for pipe flow, =hd/k 
for impingement flow 

n correlation exponent (eqs. (37), (39), and (48)); 

number of impingement holes in a single span- 
wise row in leading edge cavity 
P pressure 

Pr Prandtl number 

( P wetted perimeter 

Pt center -to-center distance between turbulators 

p p center -to-center distance between pins on tri- 

angular pitch 

q + heat flux parameter = q" ( GH 0>1 ) 
q" heat flux 

R recovery factor 

Ra rotational Rayleigh number, = 0 2 (3r*T/va 

Re Reynolds number, = DG/p 
Re w modified wall Reynolds number, 

= (DG/p)(p w /p b ) 

Re$ rotational Reynolds number, = V.D 2 /v 

r tube radius 

S Rossby number, = Q D/ V m 

S total heat-transfer surface area (exposed wall 

area + pin area) 

St Stanton number 

T temperature 

t wall thickness 

w + roughness function, 

= yjl/f +2.5 ln(2e/D) + 3.75 
u' friction velocity, =\[tq7p 

V velocity 

v open volume in pin array section 

W mass flow rate 


x distance along pipe from inlet or impingement 

plate leading edge 
x x/L 

x ' =x-(l/2)(x n /L) 

x„, yn streamwise or span wise jet hole spacing 
x p , y p streamwise or spanwise pin fin spacing 
z distance between impingement jet hole and 

target surface 

a thermal diffusivity, =k/pc\ included angle in 

triangular duct 

13 expansion coefficient 

8 flow angle of attack 

f function defined in eq. (41) 

6 rib shape angle 

p dynamic viscosity 

v kinematic viscosity 

p density 

r wall axial temperature gradient 

t 0 apparent wall shear stress 

(p correlation function for impingement heat 

transfer (eq. (33)) 
fi angular velocity 

{ } as a function of 

Subscripts: 

a based on velocity arriving at the target plate 

aw adiabatic wall 

b at bulk temperature 

c cross flow 

cp constant properties 

d based on impingement hole diameter 

eq evaluated using equivalent slot width 

/ at film temperature 

m mean 

max maximum 

min minimum 

mod modified 

n nozzle 

/ initial condition 

j jet 

/ at arc length from leading edge stagnation line 

o stagnation condition 

p pin fin 

r reference 

s stationary; static 

t turbulent 

tot total 

w wall 
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x distance from duct inlet; streamwise direction 

y spanwise direction 

z direction normal to impingement target surface 

Superscripts: 

nx exponent associated with geometry in stream- 

wise direction 

ny exponent associated with geometry in spanwise 

direction 

nz exponent associated with geometry normal to 

impingement surface 


Passages of Air-Cooled Blades 

For those in the heat -transfer community not directly 
involved with the cooling of gas turbines, a brief descrip- 
tion of the general types of cooled passages and methods 
may be helpful. Typical turbine airfoils incorporate 
passages that provide cooling by a number of methods. 
Figure 1(a) shows an example of an airfoil design that 
uses radial passages for simple convective cooling. The 
passages may be circular, rectangular, or triangular in 
cross section. 

A more complex internal geometry may incorporate 
several cooling schemes. The midchord region of the 
airfoil shown in figure 1(b) is convectively cooled by 
radial flow of cooling air through generally rectangular 
passages. The cooling air generally makes a number of 
return flows in the adjacent passages before exiting the 
airfoil. The radial flow channels are shown in the figure 
as having smooth walls. Generally, however, raised ribs 
normal, or at some angle to the direction of airflow are 
provided on the pressure and suction surfaces of the 
airfoil. These ribs, often refered to as turbulaters, break 
up the boundary layer and increase heat transfer. This 
type of airfoil design is often referred to as the multipass 
cooled configuration. At the trailing-edge cooling air- 
flows through chordwise channels before exiting the air- 
foil. Often arrays of short cylindrical or eliptical pins are 
added in this region to increase structural strength as well 
as to act as turbulators to increase heat transfer. 

The leading-edge area of the turbine airfoil is generally 
more difficult to cool than other areas because of the very 
high heat-transfer rates encountered in the gas stream at 
the stagnation area. Impingement cooling (e.g., fig. 1(b)) 
is often used. In this cooling method cooling air from the 
forward flow channel flows through single or multiple 
rows of holes to impinge on the inner surface of the 
leading-edge wall. The cooling air is then routed within 
the airfoil or ejected from holes in the leading-edge walls 
to further protect the blade with a film of cooling air. 

Often the impingement cooling method is used at the 
midchord region of turbine vanes (fig. 1(c)). In this 


method cooling air exits from a large array of holes in a 
sheet metal insert in the center of the airfoil to impinge on 
most of the airfoil inner surfaces before exiting the airfoil 
through chordwise channels, arrays of pins, or film 
cooling holes at the trailing-edge region. 

Both nonrotating turbine vanes and rotating blades use 
this cooling scheme in its various forms or in combi- 
nations with other schemes. Although frequently not 
accounted for in the thermal designs, the effects of 
rotation must be understood to improve blade passage 
design. At the high rotational speeds encountered in air- 
craft turbines, the Coriolis forces generate secondary 
forces within the flow channels. In addition, free- 
convection heat transfer could become significant 
because of the buoyancy forces caused by the density 
gradients and high centrifugal forces generated. Typical 
rotational speeds and acceleration levels for advanced 
large gas turbines are of the order of 17 000 rpm and 
75 000 g’s, and for small turbines of the order of 65 000 
rpm and 400 000 g’s. Figure 1(d) illustrates a simple 
radial passage design in a rotating turbine blade. Since 
the pressure and suction-side walls of the turbine at the 





(a) Simple convective cooling. 

(b) Cooled by several methods, 
(c) Impingement cooled vane. 

(d) Rotating blade. 

Figure 1. -Cooled turbine airfoils. 
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midchord are shown parallel to the axis of rotation, the 
Coriolis forces would act on the fluid in a direction 
normal to these two walls. However, at the trailing or 
leading edge, where the suction and pressure surfaces are 
not parallel to the axis of rotation, the Coriolis forces will 
act on the fluid in a direction that is not normal to these 
walls.. The secondary flow patterns within the passages 
would be expected to differ depending on the location of 
the flow channels. 


Convection Cooling 

Simple convective cooling is currently used only in 
older turbines, or in applications where the turbine-inlet 
temperatures and heat load demands are low and where 
more complex cooling geometries are either physically 
impractical or economically unjustifiable. In some ad- 
vanced applications simple convective cooling can still be 
used in conjunction with external ceramic coatings, 
which shield and insulate the blade metal from the hot 
gas (ref. 2). Simple convective cooling may also be used 
in applications where the air is precooled below 
compressor-exit bleed temperatures by a heat exchanger 
or by the injection of water (and use of its latent heat of 
vaporization). 

In simple convective cooling, cooling air generally 
flows in the radial direction (fig. 1(a)). In certain 
instances, however, the cooling air is made to flow in the 
axial direction. 

Most experimental investigations of forced-convection 
cooling have been performed in long, circular tubes and 
under conditions of fully developed velocity profiles 
entering the heat-transfer test section. Such carefully 
established inlet conditions help contribute reproduci- 
bility in the experimental results. However, the coolant 
passages in turbine blades pose flow and geometry condi- 
tions quite different from those maintained in most 
forced-convection heat-transfer investigations: Velocity 
profiles are not fully developed and the passage cross- 
sectional shape is apt to be noncircular. In evaluating the 
analysis methods or correlations available for turbine 
cooling application, it is necessary to examine the exten- 
sive literature on fully developed convective heat trans- 
fer. Consequently, this literature will be cited and 
discussed herein. Utilizing entrance effect studies and 
other studies pertaining to coolant passage geometry in 
combination with fully developed predictions, it may be 
possible to obtain thermal design procedures relevant to 
turbine blade or vane geometries. In addition, re- 
commendations are made herein for further heat -transfer 
research to reduce the uncertainty in thermal design 
predictions. 


Circular Passages 

Simple convective heat transfer through circular pipes 
has been extensively investigated, giving both empirical 
correlations and analytical formulations (e.g., refs. 3 to 
29). 

Siegel and Sparrow (ref. 3) addressed the question of 
how heat transfer varies in turbulent flow under bound- 
ary conditions of uniform wall temperature and uniform 
wall heat flux. They found that, for a circular pipe, the 
turbulent heat-transfer mechanism in the fully developed 
and in the thermal entrance regions are insensitive to 
these two boundary conditions for Prandtl numbers 
equal to or greater than 0.7. Therefore, the boundary 
conditions most often used do not significantly alter the 
test results or correlations when air is used as the fluid. 

One of the earliest correlations for fully developed 
turbulent flow in smooth circular pipes and for constant 
physical properties is due to Dittus and Boelter (ref. 4): 

Nu = 0.023 Re0.8p r 0.4 (1) 

where fluid properties were evaluated at the bulk 
temperatures. 

In recent years, a correlation for constant physical 
properties due to Petukhov and Popov (ref. 5) has also 
been used: 


Nu = 


(//8)Re Pr 

1.07 + 12.7V/78(Pr2/3-l) 


( 2 ) 


where / is the friction factor, which may be calculated 
from the Prandtl-Karman equation: 


=2 log (ReV7)-0.8 (3) 

Because of its simplicity and long standing usage, the 
Dittus-Boelter correlation has remained a favorite with 
which other correlations are compared. The correlation 
still gives satisfactory results for circular passages and 
other passage geometries that are not too extreme (e.g., 
very small included angles), and for nearly constant fluid 
physical properties. (Such would be the case if the heat 
flux were moderate and if the ratio of gas to coolant tem- 
perature were only slightly different from 1.) In turbine 
cooling applications the geometry of the coolant passages 
are frequently noncircular, the turbine-inlet temperatures 
are of the order of 1650 K (2500° F) or higher, the heat 
fluxes of the order of 5.5 x 106 watts per square meter, 
and wall-to-coolant temperature ratios of the order of 2.5 
to 1. Under these conditions, the assumption of constant 
physical properties is no longer valid, and the predicted 
heat-transfer c oefficien ts can be in considerable error (of 
the order of yjTi,/T w . 
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Considerable work, both analytical and experimental, 
has been done in the area of variable physical properties 
or varying heating conditions (refs. 7 to 27). The experi- 
mental data covered both passage-length-averaged and 
circumferentially averaged axial heat-transfer coeffi- 
cients. Many of the cited references report research con- 
ducted for nuclear reactor applications where hydrogen 
or helium is frequently used and where the surface-to- 
bulk temperature ratios considerably exceed those found 
in turbine-blade designs. The assumption made in these 
correlations is that the gas is chemically frozen; that is, 
chemical reaction is not included. In general, the equa- 
tion correlating the experimental data take the form 

Nu = CReO-8p r 6^p'^ " (4) 


with an L/D or x/D dependence and where the Nusselt 
number is the circumferentially averaged local value at 
the axial distance x from the entrance. Either bulk or film 
properties can be used to obtain a correlation. The effect 
of the variable properties is accounted for by (T w /Tb) ~ n 
where n~ 1/2. and b~ 0.4. When « = 0, equation (4) 
reverts to the familiar form of the Dittus-Boelter equation. 
Equation (4) indicates that the heat-transfer coefficient 
decreases as the wall-to-bulk temperature ratio increases. 

McEligot, Magee, and Leppert (ref. 8) reported on the 
effect of wall-to-bulk temperature ratios on heat transfer 
and pressure drop. The experiments covered a range of 
wall-to-bulk temperature ratios from unity to 2.5, and 
Reynolds numbers based on bulk properties at pipe inlet 
conditions from 1 .5 X 10 4 to 2.33 x 105. The experimental 
data on average heat-transfer coefficient with properties 
based on bulk temperatures agreed favorably with Taylor 
(ref. 9) and other investigators and correlated with the 
following equation: 

Nu 6 = 0.021 Re£ 8 Prg- 4 ^p^ °' 5 (5) 

Local heat-transfer coefficients were correlated to 
about 13 percent by the following correlation, except for 
the first few diameters of length: 

Nu = 0.021 Re0-8Pr0.4^pf^ _0 ' 5 |"i + 0? 1 (6) 


for x/D >5, where Nu is the local Nusselt number at an 
axial distance x from the entrance. The equation is valid 
for 1.5 x 10 4 <Re<6.0x 105 (based on bulk properties), 
0<q+ <0.004, 1 < (T w /Tfj) <2.4, and 350 K (170° F) 
<7^ <360 K (190° F). The heat flux parameter q + 
defined as 


(GH oJ ) 


( 7 ) 


Comparison of the above correlation with their experi- 
mental data showed agreement to about 13 percent. 

Taylor conducted experiments using precooled helium 
and hydrogen, with surface to bulk temperature ratios of 
up to 8, and surface temperatures up to 2950 K (5300° R) 
(ref. 10). He concluded that the circumferentially aver- 
aged heat-transfer coefficients, as a function of the axial 
distance x from the inlet, can be correlated by the Dalle- 
Donne equation of reference 17: 


Nu* = 0.021 Re£ 8 Prg- 4 ^) 


[0.29 + 0.001 9(x/£>)] 


( 8 ) 


where the fluid properties for the calculation of 
Reynolds, Prandtl, and Nusselt numbers are all evaluated 
at the bulk temperature. Taylor reported that 90 percent 
of the hydrogen data correlated to within 10 percent, 
while 80 percent of the helium data correlated to within 
10 percent of the correlation line (fig. 2). Equation (8) is 
intended for heating in fully developed turbulent flow 
(T w /Tb> 1.0). The effect on heat transfer in the 



Figure 2. - Correlation of local heat -transfer coefficients using equa- 
tion (8). 
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developing region at the entrance of the passages is not 
included in this correlation. 

Petukhov (ref. 19), in an excellent review of heat 
transfer and friction in turbulent pipe flow, gave the 
following equation, which was obtained by simplifying 
an analytical solution for variable properties: 


Nu 6 _(T w V 
Nu C p ; 6 \ Tb ) 

where 



( 9 ) 


where a = 0 for cooling and 0.3 for heating. 

Sleicher and Rouse (ref. 20) examined a number of 
correlations for constant and variable property heat 
transfer in turbulent pipe flow. For variable property 
correlation, they used a constant property correlation 
developed earlier in the same study. Using a computer 
code, they determined the modifications necessary to best 
fit the experimental heat-transfer data of four other 
investigators (e.g., refs. 15 and 26). Sleicher and Rouse 
claimed that the following equation for gases gives an 
average deviation of only 4.2 percent for the data 
examined, with only 4 of 120 points exceeding 10 percent 
deviation: 


Nu 6 = 5 +0.012 Reg- 83 (Pr ft + 0.29)^y (10) 

where n= -log(T w /T[,) 1/4 + 0.3. The authors also 
compared the experimental data with equation (9). Visual 
examination of the figure given in reference 20 indicated 
that the average deviation of equation (9) from experi- 
mental data appeared to be approximately 10 percent. 

Simoneau and Hendricks (ref. 21) correlated several 
investigators’ experimental data for large wall-to-fluid 
bulk temperature ratios and for four different gases 
(hydrogen, helium, air, and carbon dioxide). They pro- 
posed that the experimental data for the four gases could 
be correlated by a simple equation if a specific correlation 
constant C is assigned to a specific gas: 

h=c( P v)i*D-™n±y /2 (id 

For air, the constant C has a value of 0.00420. Using 
this correlation and experimental data from references 7 
and 12 for air only, the authors obtained agreement to 
+ 11 and - 10 percent of the correlating line. 

Note that equation (11) contains no fluid transport 
property terms. Simoneau and Hendricks do not suggest 
that the physical heat -transfer process is independent of 


the fluid properties. They merely point out that equation 
(11) did correlate the experimental data over the wide 
range of gas temperatures and wall-to-fluid bulk tem- 
perature ratios considered. 

Black and Sparrow (ref. 22) studied heat-transfer 
characteristics of turbulent air flow in a circular tube with 
circumferentially varying wall temperature and heat flux. 
The imposed circumferential heat flux varied by as much 
as a factor of 2. Their analysis was based on the method 
of Sparrow and Lin (ref. 23) and assumed that the ther- 
mal diffusivities in the tangential and radial directions 
were equal. The predicted circumferential variations of 
Nusselt number and air temperature were shown to be 
higher than the experimental data. At a Reynolds number 
of 4.05 X 10 4 and a circumferential flux variation of ±30 
percent about a mean, one of the analyses discussed in 
reference 23 predicted a circumferential Nusselt number 
variation of about 36 percent, while the experimental 
data showed a variation of only about 8 percent. Sparrow 
and Lin speculated that a higher value of tangential 
thermal diffusivity than the radial diffusivity may be a 
more accurate description of the heat-transfer process. 
When a ratio of tangential to radial thermal diffusivity of 
10 is used, the predicted circumferential variations 
became overly damped and became less than the experi- 
mental variations. The authors concluded that these 
observations lend support to the use of a model wherein 
the tangential to radial thermal diffusivity is substantially 
greater than 1 near the wall, but essentially unity else- 
where in the flow. 

A comparison of the circumferentially averaged experi- 
mental data of reference 22 with the Dittus-Boelter corre- 
lation (fig. 3) showed that the data were about 8 percent 
lower. The comparison showed that even with a circum- 
ferential heat flux variation of 2, the Dittus-Boelter 
correlation provides a good prediction of the average 
heat-transfer coefficient. 

For pressure drop, McEligot, et al., (ref. 8) correlated 
the friction factor to within about 8 percent of the 
Karman-Nikuradse equation 1 for x/D greater than 30. 
For engineering calculations for x/D >30, the authors 
suggested the following approximate correlation: 



where f cp is the friction factor for adiabatic, constant 
property conditions determined from the Reynolds num- 
ber based on local bulk properties. 

Taylor (ref. 11) correlated friction coefficients for both 
laminar and turbulent flow of gases. (See fig. 4.) For 


'The Karman-Nikuradse equation differs from the Prandtl-Karman 
equation only in that the former uses the fanning friction factor and 
the latter uses the Darcy- Weis back friction factor. 
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Figure 3. -Comparison of Dittus-Boelter correlation with experimen- 
tal data on tubes with circumferential heat flux variations as large 
as 2.0. (From ref. 22.) 


turbulent flow he used 423 data points from seven 
publications (refs. 9, 10, 12 to 15, and 18). The data were 
for both heating and cooling, for five different gases, 
including air, for Reynolds numbers from 5.4 xlO 3 to 
1.87x 105, and for surface-to-fluid bulk temperature 
ratios from 0.35 to 7.35. Equation (13) correlated 90 
percent of the data points to within 10 percent of the 
correlation (fig. 4(a)). For laminar flow Taylor also 
reported that equation (14) correlated friction factor for 
four gases, including air, for Reynolds numbers up to 
2.2 X 10^, and for surface-to-fluid bulk temperature 
ratios between 1.0 and 4.1 (fig. 4(b)). For turbulent flow 



( 


0.0007 + 


0.0625 \ / T b \° 5 

R eftgoj V W 


(13) 


-01 


where Re w mol } = GD/iA w (T b /T w ) . For Laminar flow 
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(a) Turbulent flow; 423 points in turbulent region. 

(b) Laminar flow; 109 points in laminar region. 

Figure 4. - Correlation of average friction factor with density evaluated at bulk temperature and viscosity evaluated at surface temperature. (From 
ref. 11.) 
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L= 8 (14) 

8 Rc»v,mod 

Petukhov (ref. 19) gave the following analytic solution 
for air and hydrogen: 

&-(*)' ,15) 

where 

n= -0.6 + 5.6(Re;)-°- 38 


results were about 10 percent lower than the Dittus- 
Boelter correlation using hydraulic diameter (fig. 5). 

Eckert and Irvine (ref. 32) investigated pressure drop 
and heat transfer in an isosceles triangle cross section 
duct with an 11.46° apex angle. For turbulent flow their 
investigations showed that the friction factor was ap- 
proximately 20 percent below the circular duct predic- 
tions and agreed well with the theoretical predictions of 
Deissler and Taylor (ref. 30 and fig. 6). For heat transfer 
the average Nusselt number for the fully developed tur- 
bulent flow at 1 14 hydraulic diameters from the inlet was 
correlated by the following equation: 


and 


Nu = 0.0325 Re0.66 


(17) 



The pressure drop due to friction is defined as 


A P= 


fL ( P V2) b 

D 2 g 


(16) 


Noncircular Passages 

Noncircular passages are frequently encountered in 
turbine-blade cooling as well as in many other applica- 
tions. The prevailing practice is to use the hydraulic 
diameter (4A/(P) and apply it to circular pipe correla- 
tions. This procedure is widely used because of its sim- 
plicity. If the passage, or duct, geometry is not too 
extreme, this procedure yields acceptable results in 
pressure drop and average heat-transfer coefficient for 
fully developed turbulent flow. If the geometrical con- 
figurations become extreme, such as those of triangular 
cross section, with very small apex angles, then both the 
pressure drop and average heat -transfer coefficient could 
deviate unacceptably from the circular duct correlations. 

The literature is fairly extensive on flow through non- 
circular ducts (e.g., refs. 30 to 53). In general, the investi- 
gators indicated that the heat transfer in noncircular 
ducts would be less than that in circular ducts under the 
same conditions. Pressure drop in noncircular ducts 
appears to be equal to or less than that in circular ducts at 
near constant property conditions, but greater than that 
in circular ducts at variable property conditions. 

Deissler and Taylor (ref. 30) presented an analysis of 
heat transfer and pressure drop for a square duct and an 
equilateral triangular duct. The analysis, which neglects 
secondary flows, predicts that for variable fluid proper- 
ties both the pressure drop and the average heat-transfer 
coefficient would fall below the values predicted for the 
circular duct. 

Emery, Neighbors, and Gessner (ref. 31) presented a 
numerical method for calculating turbulent heat transfer 
in a square duct and concluded that the heat-transfer 



Figure 5. - Comparison of heal transfer in a square duct with Dittus- 
Boelter equation for circular duct. (From ref. 31.) 



Figure 6. -Fully developed friction factors in 11.46° duct. (From ref. 
32.) Solid points calculated by method of Deissler and Taylor 
(ref. 30). 
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Properties are evaluated at bulk temperatures and 
where thermal conditions are fully established. 

Figure 7 compares the experimental heat -transfer data 
of reference 32 and equation (17) with the Dittus-Boelter 
correlation based on hydraulic diameter. The experi- 
mental data were only about one-half as large as the 
Dittus-Boelter correlation predicts. Eckert and Irving 
speculated that this difference may result from the 
turbulent heat diffusivity, on which the heat transport 
process primarily depends, being much smaller for 
extremely acute apex angles than for a round duct. 

Carlson and Irvine (ref. 33) made a systematic investi- 
gation of the friction factors of air in fully developed 
turbulent flow, using various isosceles triangles with apex 
angles ranging from 4 to 30. Their results show that the 
friction factors are generally lower than those predicted 
for the circular duct and that they vary as a function of 
the included apex angle (as shown in fig. 8). The authors, 
therefore, suggested a correlation in the form of a modi- 
fied Blasius equation 

f - era <l8) 

where C is a function of the apex angle (see fig. 9). 

Ainsworth and Jones (ref. 34) measured heat-transfer 
coefficients in circular, rectangular, and triangular ducts, 
using a short duration transient technique. The triangular 
duct had an isosceles cross section with an apex angle of 



Figure 7. -Average cross section Nusselt numbers at exit section. 
Properties evaluated at local bulk temperature; k w t/kD h -2A; 
x/Dfj = 114. (From ref. 32.) 



Figure 8. - Coefficient of friction for triangular ducts with small apex 
angles. (From ref. 33.) 



Figure 9. - Variation of the constant in the Blasius equation (eq. (18)) 
for friction as a function of the apex angle in a triangular duct. 
(From ref. 33.) 


10°. At 27 diameters from the duct inlet, where the tur- 
bulent flow was fully turbulent and at Reynolds numbers 
of 5.7 x 10 5 , the experimental data for the triangular duct 
agreed with Eckert and Irvine’s correlation (ref. 32) to 
within 8 percent (fig. 10). 

Altemani (ref. 35) studied the heat-transfer and 
pressure-drop characteristics of turbulent airflow in an 
asymmetrically heated, sharp-cornered, equilateral tri- 
angular duct. Two sides of the triangular duct were 
heated, and the third, unheated side made of low- 
conductivity material. Tests were conducted with negligi- 
ble property variations with hydraulic-diameter -based 
Reynolds numbers from 4.0 x 10 3 to 6.0 x 10 4 . The results 
(fig. 11) showed that the friction factors were 9 to 12 
percent lower than either the Blasius equation or the 
Prandtl-Karman equation. The experimental data were 
found to fit best the correlation due to Malak, Hejna, 
and Schmid (ref. 36) given in the following equation: 

^J' S = 2 log (Re ft-SK- 1 - 5 ) -0.8 (19) 

where K t = 0.936 for an equilateral triangular cross 
section. 

For Reynolds numbers between 1 X 10 4 and 5 X 106, / 
can be approximated by the following equation: 
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Figure 10. -Comparison of peripherally averaged Nusselt number for 
small-apex-angle triangular ducts. Short duration data; Reynolds 
number, 5.7 x 10 4 . (From ref. 34.) 
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REYNOLDS NUMBER, Re 

Figure 11.- Experimentally determined values ot friction factor com- 
pared with correlations. (From ref. 35.) 

/= (1.82 log Re -1.64) -2 (20) 

Figure 12 compares Altemani’s experimental heat- 
transfer data with the Dittus-Boelter correlation (eq. (1)) 
and the Petukhov-Popov equation (eq. (2)). For 
Reynolds numbers between l.OxlO 4 to 6.0 xlO 4 , the 
experimental heat -transfer data were about 15 to 22 per- 
cent lower than the Petukhov-Popov correlation and 27 
to 30 percent lower than the Dittus-Boelter correlation. 
As discussed previously, the Eckert-Irvine correlation for 
a sharp cornered, isosceles triangular cross-section duct 
with a 11.46° apex angle was about 50 percent below the 
Dittus-Boelter correlation. 

In an attempt to improve the correlation of the data 
and the Petukhov-Popov equation, Altemani used a mea- 
sured friction factor in equation (2) instead of calculating 
it from equation (20) which was intended for a circular 
duct. With the substitution, he found that the experi- 
mental data correlated with equation (2) between 1 and 
5 percent (fig. 13). He has suggested that this modifica- 
tion of the Petukhov-Popov equation could be used for 
Reynolds numbers from 1 x 10 4 to 5 X 10 6 . 

For pressure drop Hartnett, Koh, and McComas (ref. 
37) presented a comprehensive comparison of predicted 
and measured friction factors for turbulent flow through 
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Figure 12. -Measured, fully developed values of the Nusselt number 
compared with Dittus-Boelter and Petukhov correlations. (From 
ref. 35.) 



REYNOLDS NUMBER, Re 


Figure 13. -Measured, fully developed values of the Nusselt number 
compared with the Petukhov-Popov equation (eq. (2)) modified 
by Altemani. (From ref. 35.) 

rectangular ducts at ambient temperatures with no heat 
addition. Experimental data from other investigators 
(refs. 43 to 52), as well as from their own were used. They 
reported that for Reynolds numbers between 6x 103 and 
5 x 105, the circular tube correlation accurately represents 
the friction coefficient for flow through rectangular ducts 
for any and all aspect ratios (fig. 14). Examination of the 
figure shows that for the turbulent flow range, the 
agreement of most of the data is within 10 percent of the 
correlation. 

Lowdermilk, Weiland, and Livingood (ref. 38) mea- 
sured average heat-transfer coefficients and friction 
factors for air in turbulent flow at high heat fluxes for a 
square, a rectangle (aspect ratio = 5), and an equilateral 
triangle cross section. The surface temperatures ranged 
from 300 to 990 K (80° to 1320° F), and the Reynolds 
numbers ranged from 1 X 10 3 to 3.3 X 10 5 . They reported 
that the average heat-transfer coefficients could be ap- 
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Figure 14. -Comparison of experimental friction factors for rectangular ducts with circular duct correlation. (From ref. 37.) 


proximately correlated to circular duct results when the 
properties are calculated at the mean film temperature. 
Their experimental data were about 10 percent below the 
values predicted by the circular duct correlations with 
high heat fluxes in the turbulent regime (fig. 15). 
Attempts to correlate friction factors with the Karman- 
Nikuradse equation (eq. (3)) show only limited success 
(fig. 16 where the Karman-Nikuradse equation has been 
replotted as the Prandtl-Karman equation). The equation 
correlated the square duct data to within 20 percent (fig. 
16(a)). The data for rectangular and trianglar ducts were 
as much as 35 percent higher than the Karman-Nikuradse 
equation. 

Campbell and Perkins (ref. 39) investigated local fric- 
tion factors and heat-transfer coefficients for turbulent 
airflow in a vertical, rounded corner, equilateral tri- 
angular duct. The local wall-to-fluid-bulk temperature 
ratios ranged from 1.1 to 2.1, with Reynolds numbers 
ranging from 1.09xl0 4 to 3.7 xlO 4 at the inlet. Using 
bulk properties, a wall-to-bulk temperature ratio factor 
to account for the effect of variable properties, and an 
x/D h factor to account for entrance effects, they arrived 
at the following correlation: 

( T \-0.7 

-H*rosr] » 

for 6.0<x/£)/,< 123, 1.1 < T w /T b <2.\, and for proper- 
ties evaluated at the bulk temperature. The experimental 
data were correlated to 15 percent and, with all but a few 


points, to within 10 percent (fig. 17). The authors noted 
that their heat-transfer coefficient was about 10 percent 
lower than the corresponding circular tube results of 
other investigators, but that the effect of T w /T b was 
about the same. 

Campbell and Perkins also found that the ratio of 
friction factors based on both variable and constant 
properties can be correlated by the following equation for 
pressure drop, which included correction for entrance 
effects: 

f / T \-0.4 + (£V*) 0 - 67 

= (22) 
J w,cp \ 1 b / 

for 14.5 <jc/jO* < 72; and 1 . 1 < T w / T b < 2. 1 . The 
constant-property friction factor was found by using the 
modified wall Reynolds number as defined by 

Re* = 1^65 

Pw Pb 

The authors claimed that the friction factors are ex- 
tendable to x/D/j’s from 5 to 200. They concluded that 
the friction factor for a rounded corner triangular duct is 
about 20 percent higher than the corresponding circular 
duct value for a T w /T b of 2.0. 

Batista and Perkins (ref. 40) presented heat transfer 
and friction factors for turbulent airflow in a vertical 
square duct with moderate property variations, Reynolds 
numbers from 2.1 x 10 4 to 4.9 x 10 4 , maximum wall-to- 
bulk temperature ratios to 2.13, and x/D h ’s from 22 to 
155. They found that the experimental heat -transfer data 
can be correlated by equation (21) (formulated by 
Campbell and Perkins for rounded corner triangular 
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AVERAGE RATIO OF 

SURFACE SURFACE 

TEMPERATURE, TO BULK FLUID 

T w . TEMPERATURE, 

K (OR) y Tjj 

A 378 (680) 1.2 

S7 494 (890) 1.5 



REYNOLDS NUMBER, Re = (p f V b D b )/p f 

(a) Square duct. 

(b) Rectangular duct. 

(c) Triangular duct. 

Figure 15. -Correlation of heat-transfer coefficients for airflow in 
ducts with variable heat flux. Bellmouth entrance; inlet 
temperature, 297 K (535° R); properties of air evaluated at film 
temperature. (From ref. 38.) 



REYNOLDS NUMBER, Re = (p f V b D h )/|j f 

(a) Square duct. 

(b) Rectangular duct. 

(c) Triangular duct. 

Figure 16. -Correlation of friction coefficients for airflow in ducts 
with variable heat flux (Prandtl-Karman correlation (eq. (3)). 
Viscosity and density evaluated at film temperature. (From ref. 
38.) 


ducts). (See fig. 18). The properties used are based on the 
bulk temperature. Figure 19 shows a correlation of fric- 
tion factors for variable properties, normalized to con- 
stant property values and obtained at the same wall 
Reynolds numbers defined in the Symbols. Similar to 
Campbell’s observations on heat transfer for rounded 
corner triangular ducts, Batista and Perkins (ref. 40) 
found that heat transfer to rounded corner square ducts 
was also about 10 percent below the circular duct results. 
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Figure 17. — Local heat -transfer results based on bulk properties. 
6 <*/£>„ <125; \.\<T W /T b <2.\\. (From ref, 39.) 
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Figure 18. -Comparison of experimental heat transfer data of 
reference 40 with equation (21). 



tie friction factors for these geometries were about 20 
:rcent higher than those of the circular duct with the 
me heating conditions. 

From the above discussion it is clear that numerous 
iblications are available for average heat-transfer and 


pressure-drop correlations for generally fully developed 
flow downstream of the inlet in both circular and non- 
circular passages. Generally, the correlations for average 
heat transfer are about ±10 percent. Because measure- 
ment accuracies in even fundamental experiments give a 
data scatter about a correlating equation of ± 10 percent, 
the correlations for the circular tubes with fully devel- 
oped flow and uniform heating are considered adequate, 
except for extreme geometries. However, more funda- 
mental information on noncircular passages and periph- 
eral heating variations may be suggested on which to base 
more accurate prediction methods or codes for these 
conditions. Such research, however, would be of low 
priority for advanced aircraft turbines because the cool- 
ant passage length-to-diameter ratios are relatively small 
in aircraft turbines, so that most of the passage length has 
developing flow rather than fully developed flow. The 
information discussed in the next section is probably of 
more direct need for most aircraft engine applications. 

Entrance Region 

Many of the correlations given in the previous section 
are only valid for heat transfer in fully developed tur- 
bulent flow. In turbine cooling applications the coolant 
generally travels through tortuous paths before abruptly 
entering the passages that will distribute the cooling air 
within the turbine blades or vanes. Within this entrance 
region hydraulic and thermal boundary layers generally 
develop simultaneously. The heat-transfer and flow char- 
acteristics in this entrance region differ from those in a 
fully established turbulent boundary layer. In advanced 
aircraft turbines with attendant low aspect ratio blades, 
non fully developed flow could occupy a considerable 
portion of the cooled passage length. For this reason this 
section was intended to specifically emphasize the 
entrance region of the flow passage. 

Reynolds, Swearingen, and McEligot (ref. 54) pre- 
sented an analytical method for predicting heat transfer 
at the thermal entrance region of a circular duct down- 
stream of a fully developed hydrodynamic boundary 
layer. The solution, which used the method of Sparrow, 
Hallman, and Siegel (ref. 55), is based on a Reynolds- 
number-dependent velocity profile. Over a range of 
Reynolds numbers from 3.0x103 to 5.0 xlO 4 , the 
predicted values lie within 5 percent of the experimental 
data. The authors of reference 54 showed that heat trans- 
fer to the thermal entrance region of a circular duct can 
be correlated to within 5 percent for entrance length-to- 
diameter ratios greater than 2. 

Campbell and Slattery (ref. 57) presented an analytical 
approach for predicting the radial velocity and pressure 
profile downstream of the entrance of a circular tube 
with an assumed flat velocity profile. The results of their 
analysis as well as those of various investigators were 
compared with measurements of other investigators. 
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These comparisons showed that the predicted velocity 
profile development with downstream distance was 
reasonably good and that the predicted pressure-drop 
curve agreed to within 5 percent of the experimental data. 
Although the published literature (e.g., refs. 54 to 63) 
provide a better understanding of the physics of the flow 
and heat transfer at passage entrances, much of this 
research was performed under idealized conditions such 
as uniform velocities at the entrance and fully developed 
hydrodynamic boundary layers before heat addition and 
has little resemblance to entrances in turbine blades. A 
study by Boelter, Young, and Inerson (ref. 62) contains 
information that is probably of more immediate use to 
the turbine designer. Figure 20 (from ref. 62) shows the 
wide range of entrance geometries investigated, including 
bellmouths, gentle and sharp bends, long, short, and 
abrupt entrances. This array of entrance geometries, or 
slight variations of them, probably includes entrances of 
the types expected in turbine blades. The authors of refer- 


ence 62 present plots of experimental local and average 
heat transfer coefficients for several values of turbulent 
flow Reynolds numbers as a function of distance from 
the entrance for each of the geometries shown in figure 
20. Their results show that for some entrance geometries, 
the local heat-transfer coefficients at the entrance were as 
much as four times greater than the fully developed 
values in the downstream region. Altering the turbulence 
level ahead of most of the entrances (by the use of 
screens) had negligible effect on heat transfer for the 
conditions of their tests. The authors presented the 
following correlation for average heat-transfer coeffi- 
cients for entrance tube lengths as a function of tube 
length to diameter ratios greater than 5. 

h=h x (^+K~^j (23) 

where h is the average (circumferentially and axially) 



(b) BELLMOUTH WITH ONE SCREEN. (h) SMALL -ORIFICE ENTRANCE. CALMING SECTION. 




(c) BELLMOUTH WITH SCREEN HOLDER. (i) SHORT CALMING SECTION. (n) 45° ROUND-BEND ENTRANCE. 




Id) BELLMOUTH WITH SIX SCREENS. 


(e) RIGHT -ANGLE -EDGE ENTRANCE. 




(f) BARE, SHARP-EDGE ENTRANCE. 



(k) 45° ANGLE -BEND ENTRANCE. (o) 90° ROUND-BEND ENTRANCE. 



(I) 90° ANGLE-BEND ENTRANCE. Ip) 180° ROUND-BEND ENTRANCE. 


Figure 20. —Entrance conditions for which heat -transfer coefficients were determined by Boelter, et al., (ref. 62). 
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heat-transfer coefficient, h x is the downstream average 
heat-transfer coefficient where both the hydraulic and 
thermal conditions are fully established, and K is the 
tabulated value (table I) for eight of the entrance geo- 
metries shown in figure 20. No correlation was given for 
the local heat-transfer values. 

Since local heat-transfer values can be as much as four 
times greater than the fully developed downstream 
values, the effect on turbine metal temperature can be 
large. Further examination of the data of Boelter, 
Young, and Iverson (ref. 62) or additional data is needed 
to provide correlating equations for predicting local heat 
transfer coefficients as a function of distance from the 
inlet turbine blade and vane geometries. 

Rotational Effects 

At the high rotational speeds encountered in aircraft 
turbines, the Coriolis forces generate secondary flows 
within the flow channels. In addition, free-convection 
heat transfer becomes important because of the buoyancy 
forces caused by the density gradients and by the high 
centrifugal forces generated. 

Several investigators have studied the phenomenon of 
combined free- and forced-convection heat transfer due 
to rotation. Some have experimented with liquids such as 
water as the working fluid. This review, however, is 
limited to investigations with air. Several investigators 
have approached the problem analytically. The equations 
generated in these analyses are best solved through the 
use of computer codes. Skiadaressis and Spalding (ref. 
64) analyzed rectangular ducts, of aspect ratios 
(height/base) of 2: 1 , 1:1, and 1 :2, rotating around an axis 

TABLE I. - CONSTANTS K FOR PREDICTIONS OF 
AVERAGE HEAT TRANSFER COEFFICIENTS FOR 
PASSAGES WITH GEOMETRIC ENTRANCES IN 
BOELTER, ET AL. (REF. 62) 


Entrance conditions 

K 

Bellmouth (experimental) 

0.7 

Bellmouth with one screen: 
Ref. 62, Eq. (A4) 

Ref. 62, Eq. (A12) 

Ref. 62, Eq. (A8) 
Experimental 

0.144 Re 0 25 
1. 1 

0.128 Re 0 - 25 
1.2 

Short calming section with sharp-edge entrance 
(experimental) 

-3 

Long calming section with sharp-edge entrance: 
Ref. 62, Eq. (A15) 

Experimental 

0.067 Re 0 - 25 
1.4 

45° angle-bend entrance (experimental) 

-5 

90° angle-bend entrance (experimental) 

-7 

2.5-cm (1-in.) square-edge-orifice entrance 
(experimental) 

-16 

3.8-cm (l‘/i-in.) square-edge-orifice entrance 
(experimental) 

-7 


perpendicular to the ducts with radial coolant outflow 
and uniform heat flux through the duct walls. The 
authors indicated that for turbulent flow the heat- 
transfer and friction coefficients were increased by 
rotation and that heat transfer and flow were only weakly 
dependent on the aspect ratio. The increase in heat trans- 
fer between a passage with a 1:2 versus a 2:1 aspect ratio 
was about 10 percent. The analytical prediction showed a 
trend that agreed with the experimental data of Moore 
(ref. 69). Both references showed that on the pressure- 
side-wall (leeward side of rotation) heat transfer 
increased with rotation, whereas on the suction side wall 
heat transfer decreased with rotation. The end walls 
showed an increase with rotation. Moore’s data showed 
these changes in heat transfer to be approximately 10 
percent. 

Mori, Fukada, and Nakayama (ref. 66) provided both 
analytical and experimental data on a straight circular 
pipe rotating around an axis perpendicular to its own 
with radial coolant outflow. The analytical results pre- 
sented agreed reasonably well with the experimental data. 
The authors concluded that, in contrast to laminar flow, 
the local Nusselt number for fully developed turbulent 
flow had little variation in the circumferential direction, 
and that the increase in the Nusselt number over the non- 
rotating case was a little over 10 percent. 

Lokai and Limanski (ref. 70) presented experimental 
heat-transfer data in radial channels of rotating turbine 
blades with radial coolant outflow. The rotational speed 
varied from 0 to 5000 rpm. The experimental data were 
correlated on the basis of similarity parameters. The data 
spread, however, is large. The correlation does not in- 
clude the influence of a rotational Rayleigh number. 

Ito and Nanbu (ref. 68) obtained extensive experi- 
mental data on velocity profiles and friction factors for 
fully developed laminar and turbulent flows in a circular 
pipe rotating about an axis perpendicular to its own with 
radial coolant outflow. They developed an empirical 
equation for friction factors. For turbulent flow the 
following equation for the ratio of friction factors with 
and without rotation (f/f s ) was shown to agree well 
(within 3 percent) with experimental data for values of a 
parameter K t between 1 and 500: 

~ = 0.942 + 0.058 A?- 282 (24) 

Js 

where 

Re^ _ Q 2 Z)4/ 1/2 
K, ~ Re ~ V m D I v 

and where 0 is the angular velocity, V m the mean axial 
velocity, D the pipe diameter, and v the kinematic vis- 
cosity. For higher K t values, the following equation may 
be used: 
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=0.924 K } /20 (25) 

Js 

Skiadaressis and Spalding (ref. 64), in an analysis for a 
rotating circular tube with radial coolant outflow, ob- 
tained good agreement of their predicted velocity profiles 
and the experimental data of Ito and Nanbu (ref. 68). 
The analytical predictions of friction coefficients by 
Mori, Fukada, and Nakayama (ref. 66), however, 
showed better agreement with the experimental data of 
Trefethen (ref. 72). 

Morris (ref. 67) analyzed and obtained experimental 
data on the rotational effect of a circular pipe rotating 
either parallel or perpendicular to the rotating axis. 
Experimental data were obtained for rotational speeds of 
up to 2000 rpm (approximately 1400 g’s maximum at the 
midpoint of the pipe). Morris stated that, for a pipe 
rotating parallel to its axis, the Coriolis effects are more 
dominant in the entrance region and the rotational 
buoyancy effect is more important in the downstream 
developed flow region. Both these individual effects tend 
to improve heat transfer. For a pipe rotating perpendicu- 
lar to its axis, the secondary flows generated by the 
Coriolis force generally help the heat-transfer process. 
However, the buoyancy effect on heat transfer may either 
be positive or negative, depending on the direction of 
flow. For radially inward flow the effect is positive; for 
radially outward flow the buoyancy force tends to offset 
and eventually reverse the increases in heat transfer 
brought about by the Coriolis interaction. 

Morris and Ayhan (ref. 76) particularly stressed the 
possible overprediction in heat transfer which may result 
from the use of theoretically and experimentally derived 
equations that do not include the effect of buoyancy. 
Figure 21 shows that the average Nusselt number, relative 
to the nonrotating values, can very often be less than 
unity for the case of radially outward flow. Figure 21, 
which also compares the predicted results of other inves- 
tigators (refs. 65, 66, 70, and 71), shows that the 
predicted values differ considerably among the 
investigators. 

Morris and Ayhan (ref. 76) correlated experimental 
data such as shown in figure 21 for two test sections A 
and B. The range of test conditions is shown in table II. 
The resulting correlating equation for radially outward 
flow is given as follows: 

/ Ra,,\ -0.186 

Nu = 0.022 Re0 8 f J S°- 33 (26) 

where Re is the coolant flow Reynolds number, Ra is the 
rotational Rayleigh number (based on wall-bulk tempera- 
ture difference), and S is the Rossby number. 

A plot of Morris’ and Ayhan’s experimental data 
against this correlation for radially outward flow is 
shown in figure 22, where the mean scatter band for the 
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Figure 21. - Experimentally determined average Nusselt numbers (nor- 
malized to nonrotating values) for a tube rotating about its 
longitudinal axis for radially outward flow. Rotational speed, 
1000 rpm. (From ref. 76.) 


TABLE II. - NOMINAL RANGE OF VARIABLES COVERED IN 
EXPERIMENTS REPORTED IN REFERENCE 76 


Test section 

A 

B 

Length to diameter ratio 

20.62 

10.00 

Eccentricity parameter 

126.19 

65.60 

Nominal Reynolds number 

5000-15 000 

5000-15 000 

Nominal Rossby number 

Ixl0~ 2 -6xl0 -2 

4xl0- 2 -2.3xl0- 

Nominal rotational 
Rayleigh number 

2x 10' 1 -6x10 s 

5 x 10 5 - 9x 10 6 



Figure 22. -Comparison of equation (26) with experimental data for 
radially outward flow. (From ref. 76.) 
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correlation is ±15 percent. They noted that extra- 
polations of the correlating equation to the conditions of 
industrial and aircraft gas turbine engines .with radial 
coolant outflow show reductions of Nusselt numbers of 
35 and 46 percent, respectively, compared with those 
predicted for nonrotating passages. They further specu- 
late that this may explain why actual blade temperatures 
are higher than those predicted from nonrotating cascade 
tests. 

Recently Morris (ref. 77) presented a comprehensive 
treatise on heat transfer in rotating channels which in- 
cluded many of his earlier works. Figure 23 (from ref. 77) 
showed a plot similar to that of figure 21 , but for radially 
inward flow. The figure showed that for radially inward 
flow heat transfer is enhanced by rotation and increases 
with increasing Rayleigh numbers. Predicted values from 
other investigations (refs. 65, 66, 70, and 71) were com- 
pared as in figure 21. These investigations did not indi- 
cate a difference in the predicted heat-transfer values for 
inward and outward flow. 

Reference 77 gave the following equation for radial 
inward flow, which is analogous to the radially outward 
flow equation (26): 

Nu = 0.036 Re 0 - 8 (|!f) °‘ U2 s~o.0S3 (27) 



Figure 23. -Experimentally determined mean Nusselt numbers (nor- 
malized to nonrotating values) for tube rotating about its 
longitudinal axis for radially inward flow. Rotational speed, 
1000 rpm. (From ref. 77.) 


Figure 24 (from ref. 77) shows experimental data com- 
pared with equation (27). The stated mean scatter for this 
correlation is ±13 percent. 

Although some excellent analytical work has been done 
and experimental data obtained at low centrifugal and 
buoyancy conditions, the effect of centrifugal and buoy- 
ancy forces at high heat flux and high rotational speeds 
still requires verification. These buoyancy forces can 
negate or augment the heat -transfer process in the cooled 
passages (ref. 77). Typical rotative speeds and accelera- 
tion levels for advanced, large turbine engines are of the 
order of 17 000 rpm and 75 000 g’s and for small gas tur- 
bines, 65 000 rpm and 400 000 g’s; whereas the experi- 
mental data obtained were for rotative speeds from 0 to 
2000 rpm and a maximum acceleration level of 1400 g’s. 
Although empirical correlations have been obtained for 
the average heat-transfer coefficient in a tube for both 
radial coolant outflow and inflow, further research is 
needed to expand the analytical predictive capabilities, to 
extend the experiment range to higher Rayleigh and 
Reynolds numbers, and to develop reliable, simplified 
correlations suitable for design purposes. In addition, 
since many coolant passage designs incorporate tabula- 
tors (discussed later) to increase the heat-transfer 
coefficients, the interaction of the effects of turbulators 
and rotation, particularly near the passage entrance and 
exit region, needs to be investigated. Research is also 
needed to determine the effect on heat transfer when the 
coolant passage is not radial to the axis of rotation, such 
as occurs at the leading and trailing edges of the airfoil. 

Because local metal temperatures are important to 
turbine-blade life, research is also needed to obtain pre- 
diction methods and correlations for local heat-transfer 



Figure 24. — Comparison of equation (27) with experimental data for 
radially inward flow. 
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coefficients along and around the rotating coolant pas- 
sages, as well as for coolant pressure drops. Accurate 
experiments are needed to help develop the prediction 
methods and correlations. 


Impingement 

Experience has shown that, in advanced aircraft tur- 
bines, with attendant high gas temperatures and high 
pressures, simple convective cooling, if not combined 
with other forms of cooling, is not adequate to cool the 
blade material to an acceptable level. Impingement cool- 
ing is often used to cool the leading edge as well as the 
midchord areas. 

Leading Edge 

The leading-edge area in the airfoil is more difficult to 
cool than other areas because of the very high heat- 
transfer rates encountered in the gas stream at the stagna- 
tion area. Impingement cooling is often used to cool the 
leading edge. Figure 25 is a schematic sketch of an im- 
pingement cooled leading edge. Many investigators have 
studied single-row jet impingement on the inside of semi- 
cylindrical surfaces to model the leading edge of turbine 
blades (e.g., refs. 78 to 90). 

Chupp, Helm, and McFadden (ref. 78) obtained 
experimental results on geometrical configurations and 
the range of variables that are representative in aircraft 
turbines. They reported that for Reynolds numbers 
bewteen 3.0 x 103 and 1.5 x 10 4 , the following correlation 
for stagnation line heat transfer fits the experimental data 
to within 9 percent: 

Nu 0 = 0.44 Re 0 - 7 

xe*p[-0.85(!)(;Q (£)""] (28) 


(f)“ 



Figure. 25. - Schematic of impingement-cooled leading edge. 


where d is the jet hole diameter and D is the leading edge 
inside diameter. For average heat-transfer coefficient in 
the leading-edge area, the following correlation was re- 
ported to fit experimental data to within 8.7 percent: 

NU-0.63R e«.7(^) 0S (^) 06 

x-[— (oarer] <29> 

The local heat-transfer coefficients away from the stag- 
nation line can be obtained from figure 26 (from ref. 78). 
It gives the normalized Nusselt number as a function of 
distance from the stagnation line, using the distance from 
jet nozzle to impingement target as a parameter. 

References 79 to 81 represent some of the works by 
Metzger on leading-edge impingement cooling. Metzger, 
Yamashita, and Jenkins (ref. 80) reported that the maxi- 
mum heat-transfer coefficient at the leading-edge area 
(obtained at optimum nozzle-to-target distance defined in 
ref. 80) is given as 

St = Re0.27(_j =0.355 (30) 


The heat-transfer coefficient is an average over an area 
l/b, where l is the arc length distance along the semi- 
cylinder from the stagnation line and b is the equivalent 
slot width. The authors reported that the above correla- 
tion holds for Reynolds numbers (based on equivalent 
slot width) from 1150 to 6300, and for spanwise hole 
spacing to diameter ratios c/d from 1.67 to 6.67. This 
correlation does not include the effect of any variation in 
the diameter of the target cylinder. 

Damerow, et al., (ref. 82) determined the flow charac- 
teristic for rows of impingement holes in terms of dis- 
charge coefficients as a function of the jet-to-target 



DISTANCE FROM STAGNATION LINE, l/d 

Figure 26. - Normalized Nusselt number as function of distance from 
stagnation line. Impingement hole diameter, 0.041 cm (0.016 in.); 
center-to-center distance between holes, 0.328 cm (0.125 in.); 
leading-edge inside diameter, 0.239 cm (0.094 in.). 
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spacing to the jet diameter ratio z/d , and coolant passage 
Mach number M. The equation which correlated the 
experimental data within about 15 percent was: 

( 7 \ 0.045 

| J M 0 057 (31) 

where M is the coolant Mach number. 

A total-pressure loss coefficient ,Sf tot was also deter- 
mined; the correlation of the data to within about 3 per- 
cent is a function of coolant Mach number and not influ- 
enced by the flow variations as the coolant is discharged 
through the impingement holes or by the impingement 
hole spacings. The correlation given is 

tf tot = 0.97 + 0.465 M 2 (32) 


and 


pv 2 


APtot — *tot 2 


The coefficients obtained were used in a compressible 
flow network computer code to predict the flow distribu- 
tion in a simulated complex turbine-blade cooling design. 
The authors of reference 82 stated that the predictions 
compared favorably with experimental data (within 6 to 
17 percent, depending on location). 

In spite of the excellent work already done, there is still 
considerable variance among the correlations available in 
literature, as shown in figure 27 (from ref. 83). In this 
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Figure 27. — Comparison of heat -transfer correlations for impinge- 
ment into cavity using the equivalent slot width b = ird 2 /4c for 
both Nusselt and Reynolds numbers. (From ref. 83.) 


figure both the Reynolds and Nusselt numbers have been 
redefined using an equivalent slot width (6 = irrf 2 /4c) as 
the characteristic length. Table III (also from ref. 83) 
shows the restated heat-transfer correlations using the 
equivalent slot width. There are also uncertainties in 
certain areas such as at the stagnation point that need to 
be investigated. Very little has been published on the 
effect of the shape of the leading-edge cavity (other than 
cylindrical) on heat-transfer coefficient. Only limited 
data are available on the use of multiple rows of holes. 
Additional data on local heat-transfer coefficients away 
from the stagnation line would be helpful. Additionally, 
the range of variables, including the range of impinge- 
ment hole spacings and Reynolds numbers, needs to be 
extended. 

More work is also needed for cases where the spent air 
is removed through bleed holes. In advanced aircraft tur- 
bines with high gas inlet temperatures and pressures, 
impingement cooling is often combined with film cool- 
ing. The spent air, instead of exiting at the open end of 
the channel, exits through film coolant holes in the blade 
leading edge to form a film of cooling air on the external 
surface for added protection to the blade metal. 

Midchord 

Impingement cooling at the midchord region usually 
uses rows of round holes. Here, the curvature of the 
blades is small so that jet impingement on flat plate 
correlations are generally used. Figure 28 is a schematic 
sketch of an impingement cooled flat plate. Fairly exten- 
sive literature is available for flat plate impingement 
cooling. References 91 to 110 are some representative 
publications. 

Reference 91 is among the more frequently cited 
publications for jet impingement on flat plate. In these 
tests Kercher and Labakoff used a series of square arrays 
of round jets. The spent air is allowed to flow across the 
plate, creating a cross flow condition that is similar to 
those encountered in actual turbine-blade designs. The 
experimental data were correlated with the following 
relationship: 


/ t \ 0.091 

Nu = <^2 Re^Prl/3 f ± j (33) 

where the Nusselt and Reynolds numbers are based on 
the jet impingement hole diameter, <p\ and the power m 
are functions of impingement hole geometry and 
Reynolds number, and <p 2 is an attenuation factor to 
account for crossflow effects. Colladay (ref. 105) curve- 
fitted the experimental data and presented them in the 
following form: 

m = a ‘( , jy +b ‘( x z) +c ' <34 > 
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TABLE III. -RESTATED HEAT TRANSFER CORRELATIONS FOR IMPINGEMENT INTO A CAVITY-USING EQUIVALENT SLOT WIDTH 
b = rd 2 /4c FOR BOTH NUSSELT AND REYNOLDS NUMBERS (FROM REF. 83) 

lD h = 2bc(n- \)/[b + c(n- 1)]; D m = AUt\ 


Correlation 

Refer- 

Ratio of cylinder 

Ratio of nozzle 

Ratio of nozzle 

Ratio of nozzle 

Ratio of cooled 

Reynolds number. 


ence 

to nozzle 

spacing to 

to target sep- 

to target sep- 

surface target 




diameter. 

nozzle diameter, 

aration distance 

aration distance 

length to equiv- 




D/d 

c/d 

to nozzle 

to equivalent 

alent slot 






diameter, 

slot width, 

width. 






z/d 

z/b 

l/b 



Nui,,ki = 0.87 Pr J Ke i>,.q 

/b\» 5M / b V' 208 d 0.792 

Nu b , eq = 0.2Pr(^j (-j Re, cq 

Nu b .„ = C.36^y 38 Re° 6 “ a 

/ d \ 079 ( d \ 0126 / d \ o on! / b \ 0 285 

Nu b , q = 0,1, 3^-J (cj (z) U) 

/ bWdWdV 6 

Nu beq = °.63^j (-J 

27 (0(0 (f) ]) 

Nu bitq = 0.03 Pr'^J “^C; 1 

/bV j! 0 65 /A- 0 -* z 

Nu b-=°- 26 u) R C( b -j f - b -> 7 

Nu b , q = 0.,7f^V V« for ^ >7 













Figure 28. - Schematic of impingement-cooled flat plate. 


<Pl 


-exp at (*y + b i (3f)+c 1 


(35) 


1 

<P2~ 

1 +a 3 *k 


(36) 


where = (G c /Gj) ( z/d ) and where the constants are 
given in table IV. The range of variables investigated 
were Reynolds number from 3xl0 2 to 3x104, j e t 
spacing from 3.1 to 12.5 diameters, and jet to target 
spacing from 1.0 to 4.8. Figure 29 shows a comparison of 
experimental data from reference 91 with predicted 
values obtained from equation (33). 

More recently, Florschuetz, et al., (refs. 106 and 107) 
made a thorough investigation of jet impingement on a 
flat plate using both staggered and in-line jet arrays. The 
experimental variables varied from average jet Reynolds 
numbers of 5 x 103 to 5 x 104, w ith individual jet 
Reynolds numbers as low as 2.5x103 and as high as 
7 X 10 7 , jet spacing from 5 to 15 diameters in the stream- 
wise direction and 4 to 8 diameters in the spanwise 
direction, hole diameters from 0.0635 to 0.762 cm (0.025 
to 0.300 in.), and target to hole spacing ratios from 1 to 
3. The conclusion of this investigation was that the in-line 
jet arrays yielded higher heat-transfer coefficients than 


TABLE IV. -CORRELATION 
CONSTANTS FOR 
EQUATIONS (34) TO (36)“ 


[Impingement heat transfer on flat plate] 


Coefficient 

For Reynolds 
numbers of 
300 to 3000 

For Reynolds 
numbers of 
3000 to 30 000 

O] 

-0.0015 

-0.0025 

bx 

.0428 

.0685 

Cl 

.5165 

.5070 

a 2 

0.0126 

0.0260 

b 2 

-.5106 

- .8259 

Cl 

- .2057 

.3985 

«3 

0.4215 

0.4696 

*>3 

.580 

.965 


“From ref. 105. 



Figure 29. - Comparison of experimental data (ref. 91) and predicted 
Nusselt numbers. 


the staggered configurations. At some point downstream 
of the point of injection, the degradation of impingement 
heat transfer is more than made up by the convective heat 
transfer. 

As a result of this investigation, the authors re- 
commend the following correlation for impingement heat 
transfer on a flat plate: 

Nu = fl Re'" |^1 -(S,~ ^?j"p r i/3 (37) 

where Nu is the local Nusselt number in the streamwise 
direction, but averaged in the spanwise direction, and (2, 
(B, m, and n are each functions of geometry and take the 
following form: 

ororar 

The constants C, nx, ny, and nz, are each a function of 
Ct, (B, m, and n, for both in-line and staggered arrays and 
are given in table V. 

Out of 1400 data points for the in-line jet arrays, 95 
percent fall within 11 percent of the fit line, and 99 
percent fall within 16 percent of the fit line. For the 
staggered arrays with 680 data points, 95 percent of the 
points fall within 12 percent of the fit linend 99 percent of 
the points fall within 16 percent of the fit line. 

In an alternative correlation, also from Florschuetz, 
et al., the exponents m and n are held constant. It is more 
convenient to use and is nearly as accurate as equation 
(37). 
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TABLE V. - CORRELATION CONSTANTS FOR 
EQUATION (38) a 


[Impingement heat transfer on flat plate] 



Inline pattern 


Staggered pattern 


C 

nx 

ny 

nz 

C 

nx 

ny 

nz 

Q 

1.18 

-0.944 

-0.642 

0.169 

1.87 

-0.771 

-0.999 

-0.257 

m 

.612 

.059 

.032 

-.022 

.571 

.028 

.092 

.039 

u3 

.437 

-.095 

-.219 

.275 

1.03 

-.243 

-.307 

.059 

n 

.092 

-.005 

.599 

1.04 

.442 

.098 

-.003 

.304 


a From ref. 107. 


where Nuj is the Nusselt number for the first row of 
impingement jets (without the effect of crossflow) and is 
given as follows: 



/ x \ -0.554 / v \ -0.422 
Nu, = 0.383 (“^) (&) 

x ^ ^y , ' 068 R c 0.72 7 p r i/3 (40) 

The constants for equation (39) are given as follows: 



C 

nx 

ny 

nz 

n 

Inline 

0.596 

-0.103 

-0.380 

0.803 

0.561 

Staggered 

1.07 

-.198 

-.406 

.788 

.660 


Figure 30 shows a comparison of equation (40) with the 
older correlation of references 91 and 97, for two square, 
in-line arrays, that fall within the the same range of test 
conditions. The agreement with reference 91 is good, but 
there is considerable variance with reference 97. This may 
be due to differences in the test geometry and test condi- 
tions, which were not clearly given in reference 97. 

Damerow, Murtaugh, and Burggraf (ref. 82) obtained 
experimental data for discharge coefficients for impinge- 
ment on a flat plate with crossflow. The data showed the 
effects of the variables of the tests to be insignificant, and 
the authors recommended a constant discharge coeffi- 
cient of 0.78 as the best fit to the data. The spread in the 
discharge coefficients ranged between 0.7 to 0.9. 

Florschuetz and Isoda (ref. 108) reported on the flow 
distribution and discharge coefficient effects for jet array 
impingement with initial crossflow. The authors reported 
that, for the geometries tested, the discharge coefficient 
remained essentially constant when there is no initial 
crossflow. With initial crossflow, the effect on the dis- 
charge coefficient is not significant if the ratio of initial 


Figure 30. - Comparison of equation (40) for first row of impinge- 
ment jet on flat plate with references 91 and 97 for two square, 
in-line jet arrays. (From ref. 107.) 

crossflow to total jet flow ( G c /Gj ) > 1 , and the parameter 
(y„/d) ( z/d ) >12. 

This discharge coefficient decreases as (yjd ) (z/d) de- 
creases. The results of their experimental data are curve- 
fitted and presented in algebraic form in table VI. The 
function $(G C /Gj) appearing in table VI is defined as 
follows: 

C d ={(G c /Gj) = ( Gc/G . + 0 806)0.0602 (41) 

Florschuetz and Isoda also presented a flow model to 
calculate G c /Gj. For constant discharge coefficient C d 
and for negligible friction (f= 0), G c /Gj can be expressed 
as 

G c _ 1 (1 +M)sinh Bx — M sinh(l -x') 

Gj 2C d (1 +M)cosh Bx—M cosh 5(1 — x) ' '' 

where B = yJlAQC d L/z, x=x/L, x' =x-(\/2 )(x n /L), 
and M= W c / Wj, the ratio of initial crossflow to total jet 
flow. 

When discharge coefficient C d is not constant or when 
/is not equal to zero, or both, the more general equation 
given in reference 108 must be solved numerically to 
obtain G c /Gj using an equation and boundary conditions 
also given in reference 108. 

Similar to impingement cooling of the leading-edge 
area, the midchord region can use film cooling in con- 
junction with impingement cooling. The spent air from 
impingement would exit through the film cooling holes. 
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TABLE VI. -DISCHARGE COEFFICIENTS FOR IN-LINE 
JET ARRAYS WITH CONFINED COOLANT DISCHARGE 
IN ONE-DIRECTION* 


Configuration 
x„/d, y n ld, z/d 

Equation 

(b) 

Range of G c /Gj 

(5,4,1) 

c ._ 0 85 
c d- i( G c /Gj) 

0 to 0.63 
>0.63 

(5,4,2) 

n 0.85 

1.08 nO c /Gj) 

0 to 0.83 
>0.83 

(5,4,3) 

n 0.85 

l.ll £(G C /Gj) 

0 to 0.90 
>0.90 

(5,8,1) 

0.80 

C d = - 0. 169 (G c /Gj) + 0.893 
UG C /Gj) 

0 to 0.54 
0.54 to 1.5 
>1.5 

(10.4.1) 
and 

(10.8.1) 

0.76 

Q= - 0. 128(G C /Gy) + 0.825 
t(G c /Gj) 

°0 to 0.54 
c 0.54 to 1.8 
c >1.8 

a From ref. 108. 


h HO c /Gj) defined by eq. (41). 
inferred from test data. 


Hollworth, et al., (refs. 102 and 103) investigated jet 
impingement on a flat plate with vent holes on the target 
plate. The vent holes are placed either in-line with respect 
to the impingement jets or in a staggered geometry with 
respect to the impingement holes. The experimental data 
indicated that for the staggered-hole geometry, venting 
spent air through film cooling holes consistently yield 
higher heat -transfer rates (about 20 to 30 percent higher) 
than do the same arrays with edge venting of the spent air 
to one side (fig. 31). 

In view of the growing use of film-cooling in conjunc- 
tion with impingement cooling, additional research is 
needed in coupling these two cooling methods. Other 
areas that need attention are the determination of the 
significance of the effects of airfoil curvature and real 
turbine wall-to-coolant temperature ratios on correla- 
tions currently available. Research is also needed to 
determine the two dimensional heat transfer over the 
impinged surface because of the concern over possible 
large temperature gradients in the wall. 


Heat-Transfer Enhancement 

Turbulators 

The use of artificially roughened surfaces in the form 
of repeated ribs, or turbulators (fig. 32), has gained 
interest in recent years as a means of increasing heat 
transfer in the coolant passages of airfoils. The purpose 
of these ribs is to break up the boundary layer and 
thereby reduce the resistance to heat flow. This increase 
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RATIO OF DISTANCE BETWEEN IMPINGEMENT HOLE AND TARGET 
SURFACE TO HOLE DIAMETER, z/d 

Figure 3 1 .- Comparison of impingement heat transfer with staggered 
hole arrays: film-cooling-hole venting versus edge venting. (From 
ref. 102.) 



Figure 32. -Characteristic dimensions of typical turbulater design. 
Relative rib spacing, p,/D. 


in heat transfer is, of course, also accompanied by an in- 
crease in pressure loss. Several publications have ap- 
peared in literature (refs. Ill to 118), and a few of them 
have direct application to turbine airfoil design. Webb 
(ref. 112) reported repeated ribs inside a circular tube and 
indicated that the optimum heat transfer occurs at a 
pitch-to-rib-height ratio of about 10. 

For heat transfer, Webb gives the following correlation 
for Stanton number: 
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St = 


// 2 

1 +yff/l\g{e+ } Pr0.57 -u+[e+, p t /e}\ 


and friction factor / from 

Vf = 2 - 5 in ( S ) +w ^ +) ~ 3 - 75 


(43) 


(44) 


Figures 33 and 34 show the heat-transfer and pressure- 
drop experimental data and suggested correlations taken 
from reference 112. Examination of the figures showed 
that the data are correlated to within a range of 10 to 20 
percent. Webb, et al., correlated the data in figure 33 
with the following equation: 


g=4.50(e+)°- 28 


(45) 
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Figure 35. -Comparison of short pin fin data with long pin data. 
(From ref. 123.) 


The authors further stated that by using £> eq =£> — e in the 
/ and Re calculations, the data scatter in figure 34 is 



Figure 33. -Heat-transfer correlation for repeated ribs inside circular 
tubes. For e + >25, ge + =4.50 (e + ) 0 - 28 . (From ref. 112.) 
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Figure 34. -Friction correlation (eq. (47)) for repeated ribs inside cir- 
cular tubes. (From ref. 112.) 


reduced to ±6 percent. The authors recommended a 
procedure for solving the pressure-drop and heat-transfer 
coefficients that requires using e+ as an independent flow 
variable. This procedure requires the computation of St, 
/, and Re for several values of e+ and then draw a curve 
of St and / versus Re. In the fully turbulent region 
(e+>35) and for I0<p t /e<40, the St and / can be 
calculated from the following equations (ref. 112): 


St= 7 =r— y - - , (46) 

1 + V/72[4.5(e+)°- 28 PrO.57_o.950t7,/e)O.53j 

and 

vT= 2 - 5 i ”(l)" 3 ' 75+o ' 95 (?) OM <47) 

Han, Glickman, and Rohsenow (ref. 113) studied 
repeated ribs between parallel plates and determined the 
effects of rib shape, angle of attack, and pitch to height 
ratio on friction factor and heat transfer. For heat trans- 
fer, they give the following correlation: 


10(e+/35)' 
(8/45° )J 


(48) 


where 


/ = 0 when e+ <35; / = 0.28 when e+ >35 
y = 0.5 when 5<45°;y= -0.45 when 6 >45° 

For pressure drop, the following correlation is given: 

4.9(e+/35)'” 


u+ = 


( 0/9O°)°- 35 [ lO/^/e) ] * (6/45 °) 0 - 57 


(49) 


where 
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m = -0.4 if e+ <35; m = 0 if e+ >35 

«= -0.13 if p t /e< 10; n = 0.53(5/90°) if p t /e > : 10 

When 5 and 0 are both equal to 90°, e+>35, and 
p t /e> 10, this pressure-drop correlation is equivalent to 
the one given in reference 112 (eq. (47) herein), with the 
exception that the constant is changed from 0.97 to 0.95. 

In view of the increasing usage of the repeated rib, 
additional work is required to determine the effects of rib 
orientations, geometry, noncircular passages, and pas- 
sage rotation, with particular emphasis on local heat- 
transfer coefficients. Also, because the coolant passages 
with rib turbulators generally make several radial passes 
within the arifoil, reversing flow direction from pass to 
pass, the entrance effect and the turning effect are greatly 
magnified. Unpublished data showed that heat transfer 
coefficients increase by about a factor of 2 at locations 
downstream of the tightly curved returns. Additional 
research in this area is also needed. 
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Figure 36. - Comparison of long pin fin heat transfer correlation 
(eq. (54)) with short pin fin data. (From ref. 123.) 


The Reynolds number was defined as 


Pins 

Pins are frequently used in the trailing-edge region in 
turbine airfoil design. In figure 1 the chordwise ribs at the 
trailing edge are frequently replaced by pins. In such 
applications the pin height to the pin diameter ratio is not 
large (l/d< 3). These pins cannot be considered as ex- 
tended surfaces because they often cover up more surface 
on the end walls than they provide. The primary purposes 
for using these pins are (1) to add rigidity to the trailing- 
edge surface, and (2) to act as turbulators, thereby 
increasing the heat transfer. But, as with any turbulence 
generator that increases heat transfer, they also increase 
pressure drop. 

Numerous publications have appeared on pins as ex- 
tended surfaces or as tubes in cross flow (refs. 119 to 
127). They are also topics in any standard heat-transfer 
text book (e.g., refs. 28 and 29). Theoclitus (ref. 119) 
reported on heat transfer and pressure drop on many 
configurations of long tube bundles, mainly intended for 
heat exchanger usage and are of much larger h/d ratios 
than are of interest for turbine applications. Faulkner 
(ref. 121) presented a correlation for triangular pitched 
pins, but only from available experimental data, for 
longer l/d ratios. Recently, a few studies have treated 
short pins as applied to turbine-blade design (refs. 123 to 
126). Van Fossen (ref. 123) reported on heat-transfer 
characteristics on short pins with equilateral triangular 
pitch. Such pins are representative of those used in 
turbine-blade cooling design. Van Fossen suggested the 
following correlation which agreed well with experi- 
mental data: 

Nu = 0.153 ReO-685 (50) 


Re = WW l (51) 

/* 

where D' is the characteristic length ( = 4u/S) and A is 
the average flow area through the pin array section 
( A = v/L). Properties were evaluated at the Eckert refer- 
ence temperature: 

T r = 0.5 77, + 0.28 77 + 0.22 T aw (52) 

The adiabatic wall temperature used in the above equa- 
tions was calculated as follows: 

T aw =T s +R(T tot -T s ) (53) 

where R is the recovery factor ( = \/Pr). Figure 35 com- 
pares equation (50) with experimental data from refer- 
ence 123 for short fins and with the experimental data for 
pin fins obtained by two other investigators. The heat- 
transfer data for the short fins fall considerably below 
those for the long fins. Figure 36 compares experimental 
data of Van Fossen (ref. 123) with the Faulkner’s correla- 
tion, given below (ref. 121): 

BB- = [0.023 + 4. 143 exp[ -3.094 (d/p p ) 
k f 

- 0. 8903p//) 0.5075]/Re0.2946]Re0.8 Pr 1 /3 (54) 

where 
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Re= g l G °an 

/*/ 

and 

jji _ min^ 


The Faulkner correlation gives fairly acceptable values at 
high Reynolds numbers, but varies considerably at lower 
Reynolds numbers. 

Van Fossen cautions that the good correlation ob- 
tained from equation (50) may be fortuitous because the 
correlation was obtained from only two different geome- 
tries of pin configurations and the pin arrays were limited 
to equilateral triangular arrangements. There were no 
pressure-drop data reported in reference 123. 

Simoneau and Van Fossen (ref. 124) studied a single 
heated cylinder within an array of pins from one to six 
rows. A pin geometry of l/d= 3.01 for both in-line and 
staggered arrays was used. The authors concluded that 
the addition of pin arrays upstream in an in-line pattern 
of one to five rows produced an average 50 percent in- 
crease in heat -transfer level above the base case. For the 
staggered array, the average heat transfer increased by 
21, 64, 58, 46, and 46 percent, respectively, above the 
base case for one to five rows of pins upstream. 

Metzger and Haley (ref. 125) reported experimental 
results using two sets of short pins with the following 
geometry: 


curve-fitted the experimental data and obtained the 
following correlations: For x p /d- 1.5 

Nu = 0.092 Re0.707 ( 55 ) 

For x p /d- 2.5 

Nil = 0.069 ReO.728 (56) 

where Nu and Re are the averages of the local Nusselt and 
Reynolds numbers for each individual pin row in the 
stream wise direction. Fluid properties for each local 
Nusselt number are evaluated at the arithmetic average 
between the local fluid bulk and the metal end wall tem- 
peratures, and the Reynolds number is based on maxi- 
mum velocity and pin diameter. 

Figures 37 and 38 (from ref. 125) show the experi- 
mental Nusselt numbers plotted against equations (55) 
and (56), respectively. Van Fossen replotted the data 
(from figs. 37 and 38) using Reynolds number as defined 
in equation (51). The resulting plot was well correlated by 
equation (50), as seen in figure 39. 

Additional work on short pins is needed. Equations 
(50), (55), and (56) are for specific pin array geometries. 
Heat-transfer data are needed on other pin configura- 
tions and pin arrays to formulate a more general correla- 
tion. Data on presssure drop is particularly needed since 
little or no published data exist for short pins. 


Set 

x p /d 

y P /d 

l/d 

d, 

cm 

Stream- 

wise 

rows 

1 

2.5 

2.5 

1.0 

0.508 

10 

2 

1.5 

2.5 

1.0 

0.846 

10 


A duplicate set of pin arrays for each of the above geome- 
tries, but with pins made of nonconducting material, was 
also used in the experiments. The heat-transfer coeffi- 
cient increased rapidly for the first few rows, with a 
higher peak row heat transfer in the closer streamwise 
spaced array, and decreased in the downstream rows. 
They also found that nonconducting pin heat-transfer 
results, averaged over the 10 row array, followed the con- 
ducting pin heat -transfer results closely, although with a 
slight difference in Reynolds number dependence. In all 
cases the average array heat transfer was significantly 
lower than would be predicted from the larger l/d (tube 
bank) correlations. For design purposes, the authors 
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Figure 38. -Nusselt number as function of Reynolds number for x p /d= 2.5. (From ref. 125.) 



Figure 39. — Replot of data from figures 37 and 38, using characteristic dimension D' defined under equation (51). (From ref. 125.) 



Concluding Remarks 

Although published heat-transfer and, to a lesser 
extent, pressure-drop correlations are available for most 
of the general types of coolant passages used in turbine- 
blade design, areas still remain where more data are 
needed. Most of the correlations have been generated 
from data obtained at idealized inlet conditions, at low 
temperature and pressure, with uniform wall tempera- 
ture, with uniform heating, and with fully developed 
flow. Examination of available literature showed that 
accuracies of the correlations are somewhere between 
± 10 percent and ±35 percent. 

The following summarizes the results of this review on 
the status and needs for technology on heat transfer and 
pressure drop through internal passages in cooled turbine 
blades. 

Convection Through Circular and Noncircular Passages 

Correlations for fully developed, turbulent flows with 
idealized conditions, such as bellmouthed inlets or large 
plenums, are available (accuracy ± 10 percent) to predict 
average heat transfer and pressure drop through pas- 
sages. However, in actual turbine vane or blade designs, 
tortuous inlet paths and curved passage returns are the 
rule rather than the exception. Also, the passage lengths 
are such that the major portion of the flow is in the 
developing region. Although methods are available to 
account for the effects of developing flow at the inlet, 
these are mostly for idealized conditions and straight 
inlets. Little or no published experimental data are 
available for tightly curved returns, such as those used in 
the multipass passages. The little information that is 
available showed that the effect of passage curvature is to 
increase both the pressure-drop and average heat-transfer 
coefficients. For a tightly curved return typical of that 
used in turbine-blade designs, unpublished data show 
that the heat-transfer coefficients downstream of the 
return increased by a factor of about 2. Published infor- 
mation for heat transfer on the curved walls usually 
pertains to low heat flux conditions; whereas, the suction 
and pressure surface of the turbine blade walls are subject 
to high heat fluxes from the hot gas stream. Research is 
therefore needed to develop methods and correlations to 
predict local heat transfer and pressure drop in the pas- 
sage returns. Research is also needed to include the 
effects of nonuniform peripheral heating of both circular 
and noncircular passages, especially for noncircular 
passages with sharp included angles between the walls, 
where the presence of stagnant flow zones may lower the 
local heat-transfer coefficients substantially below the 
periphery averaged values and create hot spots. 


Rotation 

Although some excellent analytical and experimental 
work has been done at low centrifugal and buoyancy 
conditions, the effect of centrifugal and buoyancy forces 
at high heat flux and high rotational speeds still requires 
verification. These buoyancy forces can negate or aug- 
ment the heat-transfer process in the cooled passages. 
While empirical correlations have been obtained for the 
average heat -transfer coefficient in a tube for both radial 
coolant outflow and radial coolant inflow, further 
research is needed to expand the analytical predictive 
capabilities and to extend the experiments to higher 
Rayleigh and Reynolds numbers expected in advanced 
gas turbine engines. 

Research is also needed on the effects of rotation on 
coolant passages with turbulators, especially in the un- 
developed flow region and where the coolant passages are 
not radial to the axis of rotation (such as occurs at the 
leading and trailing edges of the airfoil). Because local 
metal temperatures are important to turbine-blade life, 
research is needed to obtain prediction methods and 
correlations for local heat-transfer coefficients along and 
around the rotating coolant passages, as well as for cool- 
ant pressure drops. High quality experiments are needed 
to help develop the prediction methods and correlations. 

Impingement 

Leading edge. - Much research has been done on im- 
pingement heat transfer in the leading-edge area. In one 
widely used reference, the experimentally obtained one- 
dimensional, spanwise averaged heat-transfer coeffi- 
cients around the leading-edge area were stated to be 
within 8.7 percent of the given correlation equation, 
within the range of the experimental variables. Yet, the 
spread in the heat-transfer data obtained by this and 
other investigators at a given Reynolds number can differ 
by a factor of 3. Additional work is needed to obtain a 
more general correlation that would reduce the apparent 
differences among the various investigators. Only very 
limited information is available on the effect of leading- 
edge cavity shape on heat transfer, of multiple rows of 
impingement holes at the leading edge, and of bleeding 
the spent cooling air through film cooling holes. Further 
research is needed in these areas to develop prediction 
methods for local heat-transfer coefficients and metal 
temperatures in the leading -edge area. 

Midchord. - Intensive research on impingement heat 
transfer on a flat plate tunnel using both in-line and 
staggered impingement hole configurations has resulted 
in a correlation of experimental, spanwise averaged 
Nusselt numbers that agreed to within 12 percent. As part 
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of planned research, the pressure-drop data will also be 
correlated for the same geometric and flow conditions as 
those used in the heat-transfer tests. The heat-transfer 
results already published and the planned pressure-drop 
correlations should satisfy the current requirements for 
one-dimensional prediction method for the turbine-blade 
midchord region. Other areas that needed attention are 
the effects of venting the spent air through film-cooling 
holes on heat transfer, the effects of curvature of the 
airfoil, and whether the effects of high turbine wall-to- 
coolant temperature ratios are significant to heat trans- 
fer. Research is also needed to determine the two dimen- 
sional heat transfer over the impinged surface because of 
the concern over possible large thermal gradients in the 
turbine blades during engine operation. 

Enhanced Heat Transfer 

Turbulators. - A number of publications have dealt 
with the effect of repeated rib type roughness on heat 
transfer and pressure drop. The correlations for both 
heat transfer and pressure drop agreed well (estimated to 
be within 10 percent) with experimental data for the 
range of experimental variables investigated. In view of 
the increasing popularity of the repeated rib design, par- 
ticularly in multipass passage blades, work is needed to 
determine the combined effects on local heat transfer and 
pressure drop of the turbulators when coupled with en- 
trance and turning effects. Other areas where research is 
needed include the effect of noncircular passages, and 
rotation, as mentioned earlier. 

Pins. - Most available correlations were obtained by 
using data from tube bundles in cross flow for heat ex- 
changers. The length-to-diameter l/D ratios of the tubes 
(or pins) are generally greater than than those used in 
turbine-blade designs for the thin trailing-edge region. A 
recent publication on four rows of short pins with length- 
to-diameter ratios of 2 and 0.5, typical of those used in 
turbine-blade cooling applications indicated the average 
heat-transfer coefficients can be lower by as much as a 
factor of 2 when compared with other data obtained 
from larger pin l/D ratios. 

The data from this limited experiment on small l/D 
pins were correlated by an equation that used only the 
Reynolds number. Apparently, the effect of l/D in the 
range of 0.5 and 2 was not measurable. The current status 
with small l/D ratio pins is such that considerable 
research is needed to obtain a clearer understanding of 
the influence of geometric variables on both heat transfer 
and pressure drop. Some research is under way in these 
areas. 

Not mentioned above, but sorely needed, is develop- 
ment of instrumentation that will provide accurate local 


measurements of flow conditions and heat transfer in 
both nonrotating and rotating environments. More accu- 
rate benchmark data are needed at low temperature and 
pressure conditions that simulate gas-turbine conditions 
for fundamental studies that contribute to the develop- 
ment of more accurate heat-transfer and pressure-drop 
prediction models and correlations. More accurate mea- 
surements are needed at high temperature and pressure 
gas conditions to evaluate and/or verify the accuracies of 
these prediction models and correlations at near engine 
conditions. 
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